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The crystal structures and phonon spectra of orthorhombic cadmium titanate with the Pbnm 
structure and of its two possible ferroelectrically distorted phases with Pbn2i and Pb2im space 
groups were calculated from first principles within the density functional theory. The obtained 
structural parameters and frequencies of Raman- and infrared-active modes are in good agreement 
with available experimental data for the Pbnm phase. Expansion of the total energy in a Taylor 
series of two order parameters showed that the ground state of the system corresponds to the Pbn2\ 
structure into which the Pbnm phase transforms through a second-order phase transition without 
intermediate phases. A substantial discrepancy between the calculated and experimentally observed 
lattice distortions and spontaneous polarization in the polar phase was explained by quantum fluc¬ 
tuations as well as by existence of twinning and competing long-period structures. 
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I. INTRODUCTION 

The calculations of phonon spectra of ten titanates 
ATiOa {A = Ca, Sr, Ba, Ra, Cd, Zn, Mg, Ge, Sn, Pb) 
with the perovskite structure in Ref. [l| have revealed that 
orthorhombic phases of CdTiOs, ZnTiOs, and MgTiOs 
with the Pbnm structure exhibit the ferroelectric insta¬ 
bility which appears in their phonon spectra as one or 
two unstable modes with symmetries and B 2 u at the 
r point. In this work, the characteristics of these modes 
are considered and the results of first-principles calcula¬ 
tions of the crystal structure, phonon spectra, and spon¬ 
taneous polarization of the parent orthorhombic CdTiOs 
phase with Pbnm structure and its two ferroelectrically 
distorted orthorhombic modifications (with Pbn2i and 
Pb2im space groups) are compared with available ex¬ 
perimental data. 

The ferroelectric phase transition in cadmium titanate 
was discovered by Smolenskii3 in 1950. Since then, 
this phase transition was studied by dielectrio^^S and X- 
ray^i^ii methods; using Raman scatteringj^i^ IR and sub¬ 
millimeter reflectance spectroscopy^^ pyroelectric current 
measurements 4 The specific features of the phase transi¬ 
tion in CdTiOs are a smallness of the lattice distortions 
and a considerable scatter in the data on the phase tran¬ 
sition temperature (50-82 K) and spontaneous polariza¬ 
tion (0.002-0.009 C/m^). Even the data on the structure 
of the low-symmetry phase are contradictory. In partic¬ 
ular, the study of the dielectric properties^ suggests that 
the 2i/m axis becomes polar, whereas structural studies 
using synchrotron radiation suggest that it is the 2i/n 
axis that becomes polar4 In Refs. H, d, and d, the anal¬ 
ysis of the temperature dependence of the dielectric con¬ 
stant and the Raman lines intensity allowed the authors 
to suggest that upon cooling one more phase transition 
accompanied by a change in the polarization direction 
occurs in CdTiOs at about 50 K. 

The paper by Fabricius and Lopez Garcfai^ is the only 
work in which CdTiOs was studied from first principles. 


In this work, the authors calculated the electric field gra¬ 
dient at cadmium atoms for different sets of atomic coor¬ 
dinates proposed in the literature for Pbnm and Pbn2i 
structures and revealed that the latter structure trans¬ 
forms into a more stable nonpolar Pbnm phase during 
relaxation. 

In view of the discrepancies between the experimen¬ 
tal data available in the literature for cadmium titanate, 
it is expedient to perform first-principles calculations in 
order to elucidate the causes responsible for these con¬ 
tradictions and to resolve them. 


II. CALCULATION TECHNIQUE 

The calculations were performed within the density 
functional theory with the pseudopotentials and the 
plane-wave expansion of wave functions as implemented 
in the ABINIT code4i The exchange-correlation interac¬ 
tion was described in the local-density approximation 
(LDA)J4 As pseudopotentials, we used the optimized 
separable nonlocal pseudopotentials 4^ generated with 
the OPIUM code, to which the local potential was added in 
order to improve their transferabilityJ4 The parameters 
used for constructing the pseudopotentials, the results 
of their testing, and other details of calculations are de¬ 
scribed in Ref. [h. 


III. RESULTS OF THE CALCULATIONS 
A. Phonon spectrum of the Pbnm phase 

The calculated frequencies of the softest phonons at 
the r point, which are responsible for the ferroelectric 
instability in orthorhombic CaTiOa, CdTiOs, ZnTiOs, 
and MgTiOs compounds with the Pbnm structure and 
in tetragonal SrTiOs with the lA/mcm structure, are 
given in Table ID In CaTiOs and SrTiOs, all phonons are 
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TABLE I. Calculated frequencies (in cm“^) of the softest 
ferroelectric modes in four titanates with the orthorhombic 
Pbnm structure and tetragonal strontium titanate. 


Mode 

SrTiOs 

CaTiOs 

CdTiOg 

ZnTiOs 

MgTiOs 

A2v, 

55 

— 

— 

— 

— 

Eu 

39 

— 

— 

— 

— 

Bzu 

— 

82 

54 

73 

115 

B2u 

— 

97 

81i 

54* 

81 

Biu 

— 

82 

104i 

103* 

133* 

Lis 

68i 

165i 

187i 

240* 

260i 


Note: The frequencies of unstable TO phonons at the T point of 
the Brillouin zone in the cubic parent phases are shown in the 
lowest row™ 


TABLE If. Lattice parameters a, b, and c (in A) and atomic 
coordinates for CdTiOs with the Pbnm structure. 


Parameter 

This work 

Experiment 



Ref. 15 

Ref. r 

a 

5.2427 

5.3053 

5.284 

b 

5.3815 

5.4215 

5.403 

c 

7.5744 

7.6176 

7.590 

Cd. 

-0.01017 

-0.00847 

-0.00891 

Cd, 

-10.04637 

-10.03873 

-f0.03997 

Cd. 

-10.25000 

-10.25000 

-f0.25000 

TU 

-10.00000 

-fO.OOOOO 

-fO.OOOOO 

Tig 

-10.50000 

-10.50000 

-fO.50000 

Ti, 

-10.00000 

-fO.OOOOO 

-fO.OOOOO 

ou 

-10.10130 

-f0.0902 

-f0.0918 

Oly 

-10.46252 

-fO.4722 

-f0.4714 

OR 

-10.25000 

-10.25000 

+0.2500 

02, 

-10.69348 

-f0.7008 

+0.70083 

02, 

-10.30304 

-fO.2969 

+0.29660 

02, 

-10.05341 

-f0.0472 

+0.04783 


* At T = 150 K. 


stable (this corresponds to experiment), whereas in the 
other three titanates one or two unstable modes (whose 
frequencies are imaginary) arise. The strongest insta¬ 
bility in these compounds is associated with the Biu 
phonon, which can induce the Pbnm —>■ Pbn2i phase 
transition. In CdTiOa and ZnTiOa, the phonon with the 
B 2 u symmetry is also unstable. This phonon can induce 
the Pbnm —>■ Pb2im phase transition. It should be noted 
that, among the three compounds under consideration, 
the Pbnm structure can be obtained only for CdTiOa 
(ZnTiOa and MgTiOa usually crystallize in the ilmenite 
structure). This is why hereafter only the properties of 
cadmium titanate will be considered. 

The calculated lattice parameters and equilibrium 
atomic coordinates for orthorhombic CdTiOa with the 
Pbnm structure are compared with available experimen¬ 
tal data in Table El It is seen that the results of the 


calculations are in good agreement with the experimen¬ 
tal data. It should be noted that a better agreement 
is observed for the experimental data obtained at lower 
temperatures. A small systematic underestimation of the 
calculated lattice parameters is typical of the LDA ap¬ 
proximation used in this work. 

The vibrational spectrum of a crystal with the Pbnm 
space group consists of 60 modes, including 24 Raman- 
active modes (with symmetries Ag, Big, B 2 g, and B^g), 
25 IR-active optical modes (with symmetries Biu, B 2 u, 
and i? 3 „), 8 so-called silent Au optical modes, and three 
acoustic modes inactive in the optical spectra. 

The calculated frequencies of modes active in IR re¬ 
flectance and Raman scattering are compared with avail¬ 
able experimental datap^i^i^ in Table IIIIl For the Raman- 
active modes, the results of the calculations agree well 
with the experimental data obtained on ceramic sam¬ 
ples^ and single crystals;^ the typical relative deviation 
of frequencies is about 3%. A comparison of the mode 
frequencies determined from polarized Raman scattering 
spectra^ with the results of our calculations shows that 
the peaks observed at 303 and 392 cm“^ in theyy polar¬ 
ization (in the crystal setting accepted in Ref. @) is most 
likely erroneously assigned to the Big modes (in our crys¬ 
tal setting). The positions of these peaks are close to the 
positions of the peaks of the Ag modes, which can also 
be observed in the above polarization. The experimental 
peak at 307 cm“^ in the yz polarization, assigned to the 
B 2 g or Bsg modes, is most likely a “leak” of the Ag mode. 

A comparison of the calculated mode frequencies with 
the results obtained from an analysis of IR reflectance 
spectra^ appeared to be a more complex problem. A di¬ 
rect comparison of the experimental frequencies with the 
frequencies calculated for the symmetries of modes indi¬ 
cated in Ref. revealed their significant disagreement. 
In order to reliably identify the experimentally observed 
modes, the oscillator strength fi was additionally calcu¬ 
lated for each mode (the oscillator strength characterizes 
the contribution of the mode into the complex dielectric 
constant) and it was supposed that in Ref. l9| the symme¬ 
tries of the observed modes were identified incorrectly. 
This can be a result of the high twin density character¬ 
istic of CdTiOs crystals which can undoubtedly 
manifest itself in measurements on samples with a large 
area. 

The calculated frequencies and oscillator strengths for 
the CdTiOs crystal with the Pbnm structure are com¬ 
pared with the results of the IR reflectance spectra stud¬ 
ies in Fig.[TJ Under the assumption that the modes iden- 
tihed as B^u in Ref. [§ for the Pnma crystal setting cor¬ 
respond to the Biu modes in our Pbnm setting and that 
the Biu modes in Ref. correspond to the B 2 u and 
modes in our setting, the agreement between the calcu¬ 
lated and experimental data becomes more reasonable 
(Fig. E Table EH- A considerable shift in frequencies of 
three softest Biu modes and two softest B 2 u modes as 
compared to the calculated frequencies is explained by 
anharmonicity. The modes with calculated frequencies 
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TABLE III. Frequencies Vi of the optical modes active in Raman and IR reflectance spectra and oscillator strengths fi for 
IR-active modes in CdTiOs with the Pbnm structure. 


Mode 

This work 

Vi (cm ^) 
Exp. (Ref. ^ 

Exp. (Ref. 

Mode 

fi-W\ 
arb. units 

Vi 

This work 

(cm ^) 

Exp. (Ref. ^ 


96 

95 

99 

Biu 

3.87 

54 

61? 


128 

123 

125 


1.27 

104 

Ill 


195 

190 

194 


0.60 

167 

165 


295 

295 

299 


0.38 

275 

284 


414 

390 

390 


0.30 

302 

306 


449 

461 

465 


0.71 

379 

383 


512 

— 

496 


0.13 

423 

— 

Big 

140 

135,141 

144 


0.17 

445 

458 


211 

— 

— 


0.96 

513 

525 


353 

342 

346 

B2u 

5.98 

81i 

61** 


486 

— 

479* 


0.62 

90 

111** 


683 

— 

— 


0.86 

177 

177 

B2g 

117 

— 

114 


0.01 

202 

— 


280 

— 

307? 


0.37 

321 

321 


445 

— 

459* 


0.14 

345 

349 


488 

— 

509* 


0.12 

437 

— 


764 

— 

— 


1.04 

476 

483 

Big 

111 

110 

115 


0.02 

530 

525 


142 

141 

141 

Biu 

3.42 

104i 

44* * 


205 

— 

— 


0.89 

60 

96** 


356 

— 

— 


0.78 

128 

143** 


448 

— 

— 


1.24 

224 

225 


492 

— 

504 


0.62 

396 

387 


738 

— 

— 


0.02 

428 

— 






1.19 

481 

511 


* Components of a broad weakly structured line. 

** A strong deviation from the calculations is a result of anharmonicity. 


of 54 and 104 cm“^ are most likely indistinguishable in 
the experiment from the other modes, and four modes 
with the smallest oscillator strengths are not observed in 
the spectra at all. 


B. Lattice distortions accompanying the 
ferroelectric phase transitions 

As shown above, the used technique of first-principles 
calculations describes well the properties of orthorhom¬ 
bic CdTiOa with the Pbnm structure. Now, we con¬ 
sider the properties of ferroelectrically distorted phases 
of cadmium titanate. The equilibrium atomic positions 
in the distorted phases were determined by relaxing the 
Hellmann-Feynman forces in structures obtained from 
the parent nonpolar Pbnm phase by adding a perturba¬ 
tion with the Biu symmetry (for the Pbn2i phase) or i? 2 u 
symmetry (for the Pb2im phase). The structure of the 
nonpolar Pbnm phase and the character of atomic dis¬ 
placements accompanying the transitions from it to the 


low-symmetry phases are shown in Fig. [5] The transi¬ 
tion to the Pb2im phase is accompanied by out-of-phase 
displacements of the titanium and oxygen atoms along 
the y axis, whereas the cadmium atoms remain almost 
unmoved. Upon the transition to the Pbn2i phase, both 
metal atoms are slightly displaced along the z axis (by 
approximately equal distances), and the displacements 
of the titanium and oxygen atoms contain a considerable 
y-component in addition to the z-component. 

The transition to the polar phases results in the ap¬ 
pearance of two nonequivalent sets of oxygen atoms 02 
(labeled by the a and b letters in Table ITVT) . and the dis¬ 
appearance of the n plane upon transition to the Pb2im 
phase results in the appearance of two nonequivalent cad¬ 
mium atoms (Cdl, Cd2). The energy gain (per formula 
unit) associated with the distortions is equal to 6.21 meV 
for the Pbn2i phase and 1.38 meV for the Pb2im phase. 

A comparison of the calculated atomic displacements 
and the lattice strains with the results of low-temperature 
X-ray measurements^ii (Table HVl) demonstrates that the 
atomic displacements and the lattice strains in the exper- 
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FIG. 1. Comparison of (a) the calculated frequencies and 
oscillator strengths of the IR-active modes for cadmium ti- 
tanate with the Pbnm structure with (b) the corresponding 
parameters obtained from an analysis of the IR reflectance 
spectra. Labels near the points indicate the mode symmetry 
(Blu, B2u, Bsu)- 

iment are much smaller. For example, from the compar¬ 
ison of the calculated lattice parameters in Tables HIl and 
m it follows that a considerable increase (by 0.045 A) 
in the lattice parameter along the polar c axis should be 
observed in the Pbn2i phase, and the a lattice parameter 
rather than the b lattice parameter should be increased in 
the Pb2im phase. In the experiment, the largest sponta¬ 
neous strain below the phase transition temperature was 
observed for the b lattice parameter (elongation of the 
order of 0.002 A)^ which disagrees with the predictions 
for both ferroelectric phases. The possible causes of these 
discrepancies will be discussed below. 


C. Parameters of an effective Hamiltonian 

Torgashev et al<^ and Gorshunov et al<^ have supposed 
that below the ferroelectric phase transition temperature, 
one more phase transition accompanied by a change in 
the polarization direction occurs in cadmium titanate. 
In order to check this hypothesis, we calculated the de¬ 


TABLE IV. Lattice parameters a, b, and c (in A) and atomic 
coordinates for ferroelectrically distorted Pbn2i and Pb2im 
phases of cadmium titanate. 

Pbn2\ Pb2\m 

Parameter This work Exp. (Ref. ^ This work Exp. (Ref. Tj 


a 

5.2392 

5.2946 

5.2498 

5.281 

b 

5.3777 

5.4151 

5.3870 

5.403 

c 

7.6192 

7.6029 

7.5699 

7.583 

CdG 

-0.01101 

-0.0083 

-0.01400 

-0.01106 

Cdl„ 

-b0.04425 

+0.0407 

+0.04583 

+0.04076 

CdG 

-bO.25324 

+0.25 

+0.25000 

+0.25000 

Cd2, 

-bO.OllOl 

+0.0083 

+0.00509 

+0.00697 

Cd2j, 

-0.04425 

-0.0407 

-0.04696 

-0.04076 

Cd2, 

-TO. 75324 

0.75 

-0.25000 

-0.25000 

TR 

-bO.00080 

+0.004 

+0.00459 

+0.00190 

Tiy 

-bO.49548 

+0.493 

+0.50699 

+0.5045 

Ti, 

-b0.00334 

+0.004 

+0.00080 

+0.00214 

Olflx 

-to. 10277 

+0.091 

+0.10313 

+0.0925 

Oluj, 

-b0.46191 

+0.473 

+0.45559 

+0.4759 

OlOz 

-b0.24044 

+0.241 

+0.25000 

+0.25000 

Olbx 

-0.10277 

-0.091 

-0.09897 

-0.0911 

Olby 

-0.46191 

-0.473 

+0.53144 

+0.5338 

Olbz 

-to. 74044 

+0.741 

-0.25000 

-0.25000 

02ax 

-bO.68779 

+0.723 

+0.69372 

+0.6999 

02ay 

-b0.31805 

+0.308 

+0.29494 

+0.2951 

02az 

-b0.04267 

+0.047 

+0.05436 

+0.0501 

026, 

-bO.29839 

+0.323 

+0.30806 

+0.2989 

02by 

-b0.71471 

+0.710 

+0.68835 

+0.7017 

02bz 

-0.06583 

-0.045 

-0.05264 

-0.0457 


* Correct signs of atomic displacements are recovered. 

pendence of the total energy of the crystal on the lattice 
distortion with the and i? 2 ii symmetries. 

As is known, the normal modes form a complete basis 
set using which any combination of atomic displacements 
can be expanded. We expanded the lattice distortions in 
the Pbn2i and Pb2im phases determined in Sec. ME 
using the basis of eigenvectors of the dynamic matrix of 
CdTiOs with the Pbnm structure. For the Pbn2i phase, 
this expansion contains the contributions of seven modes 
with the Biu symmetry (the acoustic Biu mode that de¬ 
scribes a uniform shift of the unit cell was excluded) and 
seven fully symmetric Ag modes (Fig. [3]). For the Pb2im 
phase, the expansion contains the contributions of nine 
B 2 u modes (also without the acoustic mode) and seven 
Ag modes. In both cases, the contributions of the lowest- 
frequency modes with Bi^ and i? 2 ii symmetries to the 
total energy of distortions amount to 92.8 and 95.4%, 
respectively. Therefore, when constructing an effective 
Hamiltonian, we can restrict ourselves to the expansion 
in powers of the amplitudes of these two modes. 

The total energy of the crystal was expanded in a Tay¬ 
lor series in powers of the distortion amplitudes ^ and 
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FIG. 2. (a) Projection of the structure of the orthorhombic CdTiOa phase with the Pbnm structure onto the he plane and (b, 
c) the character of atomic displacements accompanying the ferroelectric transitions to the (b) Pbn2i and (c) Pb2im phases. 



V, cm"' 

FIG. 3. Relative contributions of different modes with Biu 
and i? 2 u symmetries to the ferroelectric distortions accom¬ 
panying the transitions to the Pb2im and Pbn2i phases of 
cadmium titanate. 


r] of the lowest-frequency normal modes with Biu and 
B 2 u symmetries. In this expansion, the lattice param¬ 
eters were assumed to be fixed and equal to the lattice 
parameters of the Pbnm phase. The resulting expansion 
had the form 

rj) = Etot{0, 0) -b -b ci^'* -b 62 ? 7 ^ + C2?7'‘ + 

( 1 ) 

with the coefficients bi = —0.4075 Ha (Hartree), 62 = 
-0.2614 Ha, ci = 183.49 Ha, C 2 = 249.58 Ha, and 
d = 457.4 Ha. It turned out that the total energy is 
adequately described by the fourth-order polynomials of 
two order parameters and that there is no need to in¬ 
clude the sixth-order invariants into the expansion. This 
means that the phase transition in CdTiOa is far from 
the tricritical point (the closeness to which was supposed 
in Ref. [^. Furthermore, since we have d > 2 ^CiC 2 , the 
minima of the total energy correspond to the order pa¬ 
rameters ±(^, 0 ) and ±( 0 , r]) and are separated from each 
other by the energy barriers. This means that the for¬ 


mation of the monoclinic phase with Pb space group 
and tilted polarization vector is energetically unfavor¬ 
able. The calculation using the formula (1) shows that 
the energy minima are observed at the mode amplitudes 
of = 0.03332, rj = 0) and (^ = 0, 77 = 0.02288) and 
that the energy gains associated with these lattice distor¬ 
tions are 6.16 and 1.86 meV, respectively. They are close 
to the corresponding energy gains of the true distortions 
(see Sec. Him. 

Since only two modes with the lowest frequencies were 
included in expansion ( 1 ), it was necessary to check that 
the energetically most favorable Pbn2i phase remains 
stable with respect to small distortions with the B 2 U sym¬ 
metry when taking into account all modes involved in the 
distortion. For this purpose, we calculated the phonon 
spectrum of CdTiOa in the Pbn2i phase. In this spec¬ 
trum, the ferroelectric mode polarized along the y axis 
had the lowest frequency of 83 cm“^. The positive val¬ 
ues of all mode frequencies confirm that the polar Pbn2i 
phase is the ground-state structure of CdTiOa. 

Although the above analysis of the polar Pbn2i phase 
stability is for T = 0 , it is unlikely that a change of the 
Cl, C 2 , and d coefficients in the thermodynamic potential 
with increasing temperature can result in the violation 
of the d > 2 y/ciC 2 condition and the transition into the 
monoclinic phase. Therefore, the possibility of succes¬ 
sive ferroelectric phase transitions in CdTiOa with vary¬ 
ing temperature can be most likely ruled out. Possibly, 
the specific features observed in the dielectric constant of 
CdTiOg near 50 are associated with the existence 

of side minima of the total energy that are separated from 
the main minima by low potential barriers. 


D. Spontaneous polarization 

The spontaneous polarization Pg in the orthorhombic 
Pbn2i and Pb2im phases of cadmium titanate was cal¬ 
culated by the Berry’s phase method.— The calculated 
values of Pg corresponding to the lattice distortions de¬ 
termined in Sec. IHI 51 for the Pbn2i and Pb2im phases 
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are equal to 0.29 and 0.21 C/m^, respectively. The spon¬ 
taneous polarization Pg corresponding to the minima in 
expansion (1) is equal to 0.21 C/m^ for the Pbn2i phase 
and 0.16 C/m^ for the Pb2im phase. It is worth not¬ 
ing that both calculated values of Pg considerably ex¬ 
ceed the experimentally obtained polarizations of 0.002- 
0.009 C/m^ (Refs, i and i). 


IV. DISCUSSION 

As shown in Ref. [l|, the ferroelectric instability is char¬ 
acteristic of cubic phases of all ten titanates with the 
perovskite structure studied in the cited work. In five of 
these compounds, which undergo structural phase tran¬ 
sitions to the Pbnm or lA/mcm phases, the ferroelectric 
instability is weakened and can be observed only in three 
compounds with the Pbnm structure: CdTiOa, ZnTiOa, 
and MgTiOa (Table U). A comparison of the frequencies 
of the softest ferroelectric modes at the T point in the 
cubic parent phases and in the Pbnm or lA/mcm phases 
shows that upon the structural phase transition, the fer¬ 
roelectric instability is retained in crystals in which it is 
strongest in the cubic phase. Therefore, it is not surpris¬ 
ing that the temperature of the ferroelectric phase tran¬ 
sition in CdTiOa is considerably lower than in BaTiOs 
or PbTiOa. 

Studies of SrTiOs have long established that despite 
the existence of the ferroelectric instability in this com¬ 
pound, the corresponding phase transition does not occur 
with decreasing temperature. It is believed that such a 
behavior is due to quantum fluctuations, i.e., zero-point 
vibrations of atomsJ^ One more factor responsible for 
this behavior can be the structural phase transition that 
occurs in this compound, which, as shown above, weak¬ 
ens the ferroelectric instability. 

It is possible that the strong influence of quantum fluc¬ 
tuations on the physical properties of crystals should also 
be observed in CdTiOa. As shown in Ref. [I^, quan¬ 
tum effects most strongly affect the modes with a small 
“reduced mass,” in particular, the ferroelectric mode to 
which light oxygen atoms make a significant contribution. 
This explains why the substitution of the isotope by 
the heavier isotope in strontium titanate results in a 
real ferroelectric phase transition)^ 

In order to take into account quantum fluctuations, ex¬ 
pression (I) for the “potential” energy should be comple¬ 
mented by the kinetic energy of nuclei. When performing 
first-principles calculations, this energy is ignored in or¬ 
der to calculate correctly the forces acting on atoms (the 
Born-Oppenheimer approximation). 

Although the approach used below requires a more rig - 
orous justification, it is possible to try (following Ref.fl^ 
to take into account quantum fluctuations at the level of 
the mode motion by adding the kinetic energy of the 
mode described by the (—?i^/2M*)V^ operator, where 
M* is some reduced mass, to the potential energy (1). 
In order to determine the quantitative criterion of the 


TABLE V. Frequencies of unstable ferroelectric modes in the 
high-symmetry phases, energy gains resulting from the lat¬ 
tice distortion™ and hu/Eo ratios for three ferroelectric phase 
transitions in SrTiOa and CdTiOa. 


Phase transition 

V (cm ^) 

Eq (meV) 

hvjEo 

SrTiOa, Pm3m —>■ R3m 

68i 

0.75 

11.2 

CdTiOs, Pbnm —> Pb2\m 

81i 

1.38 

7.28 

CdTiOs, Pbnm —>■ Pbn2i 

1041 

6.21 

2.08 


complete suppression of a phase transition by quantum 
fluctuations, we consider a simple one-dimensional prob¬ 
lem in which a particle of mass M* moves in a double¬ 
well potential V{x) = —ax^ -\- bx'^. Suppression of the 
phase transition by quantum fluctuations occurs at the 
moment when the kinetic energy of zero-point mode vi¬ 
brations becomes higher than the well depth of the po¬ 
tential barrier under consideration. The numerical so¬ 
lution of the S chrodi nger equation shows that this oc¬ 
curs when h/\/M*a^ > 0.428. By eliminating the un¬ 
known mass M* from this condition, the criterion can be 
rewritten in a physically clear form Tiiz/Eq > 2.419. Here, 
iz = 2a/M*/27r is the value of an imaginary mode fre¬ 

quency in the high-symmetry configuration (maximum of 
the potential barrier) determined in the classical (Born- 
Oppenheimer) approximation, and Eq = a^/Ab is the 
depth of the potential well for the V (x) potential. For the 
four-minimum potential well described by equation (1), 
the quantitative value of the criterion can be slightly dif¬ 
ferent. 

We use the obtained criterion to evaluate the degree 
of influence of quantum fluctuations on the ferroelectric 
phase transitions in CdTiOa. The frequencies of unsta¬ 
ble ferroelectric modes determined from first principles 
in the Born-Oppenheimer approximation, the energies 
Eq of the ordered phases, and the Hv/Eq ratios for the 
hypothetical ferroelectric phase transition in SrTiOa and 
two phase transitions in CdTiOa are given in Table Ivl It 
follows from the table that quantum fluctuations should 
suppress the ferroelectric phase transition in strontium 
titanate and the transition to the Pb2im phase in cad¬ 
mium titanate. As regards the phase transition to the 
Pbn2i phase in cadmium titanate, this transition is sus¬ 
ceptible to strong quantum fluctuations but is not com¬ 
pletely suppressed. Therefore, the only phase transition 
that can be associated with the experimentally observed 
ferroelectric phase transition in CdTiOs is that to the 
Pbn2i phase. 

Quantum fluctuations can also be responsible for the 
significant discrepancies in the structural positions of 
atoms and in the values of the spontaneous polarization 
Pg. In quantum-mechanical calculations of the ground 
state for multiwell potentials, the displacement corre¬ 
sponding to the most probable atomic position is always 
smaller than the displacement corresponding to the min¬ 
imum of the potential energy. Therefore, quantum flue- 
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tuations should decrease the distortions accompanying 
the phase transition and decrease the spontaneous po¬ 
larization Pg. Other possible factors responsible for the 
decrease in the experimentally obtained polarization in 
CdTiOa can be the manifestation of twinning in crys- 
tals^ d^d^ and the presence of long-period structures com¬ 
peting with the Pbnm phasei^ 

V. CONCLUSIONS 

The first-principles calculations of the structural pa¬ 
rameters and the phonon spectrum of orthorhombic 
CdTiOa enabled to refine the identification of Raman 
scattering and IR reflectance spectra. The calculated de¬ 
pendence of the total energy of cadmium titanate on the 


amplitudes of two unstable modes indicates that the fer- 
roelectrically distorted Pbn2i phase is the ground state 
of the crystal at T = 0. This phase appears to be the 
most stable with respect to quantum fluctuations, which 
are rather strong and suppress other possible lattice dis¬ 
tortions. Quantum fluctuations were found to be one 
of the main factors responsible for the discrepancy be¬ 
tween the calculated and experimentally observed values 
of spontaneous polarization and structural distortions ac¬ 
companying the phase transition. 
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